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We investigate the correlations between any number of arbitrarily far-apart regions 
of the vacuum of the free Klein-Gordon field by means of its finite duration coupling 
to an equal number of localized detectors. We show that the correlations between any 
TV such regions enable us to distill an TV-partite W state, and therefore exhibit true 
TV-fold entanglement. Furthermore, we show that for TV = 3, the correlations cannot 
be reproduced by a hybrid local-nonlocal hidden-variable model. For TV > 4 the issue 
remains open. 



In a recent paper [[J, the nature of the correlations 
between two arbitrarily far-apart regions of the ground 
state of the free Klein-Gordon field was investigated by 
means of its finite duration coupling to a pair of local- 
ized detectors. It was shown that a local hidden-variable 
model cannot account for these correlations p], , and 
that as a function of the separation between the regions, 
L, and the duration of the coupling, T, the entanglement 
decreases at a slower rate than e _ ( L / cT ) . It is, therefore, 
natural to ask whether the vacuum admits other types of 
these kinds of correlations, i.e. true many-region entan- 
glement , and full nonlocality £| . In this paper we will 
answer the first question affirmatively for any number of 
arbitrarily far-apart regions of the vacuum, TV, while to 
the second question we will only be able to provide a 
positive answer for the case TV = 3. We will follow the 
method developed in previous papers 0, ||, and begin 
by investigating in detail the three-region case. We will 
then use results thus obtained to study the correlations 
between any number of regions. 

Our method consists of the finite duration coupling of 
the field we wish to investigate to any number of local- 
ized nonentangled detectors, such that all the detectors 
remain causally disconnected from one another through- 
out the interaction. Once the interaction is over, we trace 
over the field degrees of freedom to obtain the detectors' 
reduced density matrix. The crux of the method lies in 
the fact that since the detectors are initially nonentan- 
gled, any nonlocal correlations exhibited by the detec- 
tors' final reduced density matrix must have their origin 
in corresponding vacuum correlations. This enables us to 



apply recently developed tools from the field of quantum 
information theory to study the structure of the vacuum. 

Before we begin, let us first give the definitions of true 
multi-fold entanglement and full nonlocality. A multi- 
partite mixed state is said to be truly multi-fold en- 
tangled iff it cannot be expressed as a convex sum of 
decomposable terms. In the tri-partite case this just 
means that the state cannot be written as a convex 
sum of terms of the form pi eg) pj^, where the sub- 
scripts denote any of the three subsystems. Examples 
of truly tri-fold entangled states are the GHZ [6( and 
W |3 states. Analogously, we can also distinguish be- 
tween different types of nonlocality. A multi-partite 
state is fully nonlocal if there does not exist a decom- 
posable hidden-variable dependent probability function 
that can account for the results of any set of von Neu- 
mann measurements. In the tri-partite case of a sys- 
tem composed of three parts A, B, and C, this means 
that pabc (a, b, c | A) ^ Pa (a I A) p B c (b, c | A) + 
pc (c | A) pab (a, b | A) + p B (b | A) p C A (c, a | A) . Here 
A is the hidden-variable, and Pabc (a, b, c | A) is the 
probability for a = a, b — b, and c = c. Otherwise, the 
state may admit a hybrid local-nonlocal hidden-variable 
description Svetlichny derived a Bell-like inequal- 
ity to distinguish between these two cases y| , which was 
later generalized [jj- 

Let us consider the ground state of a free Klein-Gordon 
field and three nonentangled point-like two-level detec- 
tors [ljj. The interaction Hamiltonian of the field and 
the detectors is given by 



Hj{t) = Hf(t)+Hf(t) + H<f(t) (1) 



f dt'c l {t')(e^ t 'a+ +e-' in ' t 'a-)<P(x t ,t') : 

=A,B,C^ T / 2 



where cf> (x, t) is a free Klein-Gordon field in three spatial dimensions, the af are the detectors' "ladder" operators, 
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and the f2j denote the energy gap between detector en- 
ergy levels. T is the duration of the interaction, while 
the window- functions, ti (t), govern its strength. We set 
cT << Lij, with Lij = \xi — Xj \ and the Xi being the lo- 
cations of the detectors. This ensures that the detectors 
remain causally disconnected throughout the interaction. 
The evolution operator therefore factors to a product. In 



the Dirac interaction representation, employing "natu- 
ral" units {h = c = 1), U = Yl l=A . B . C Te~ lS dtH]{t \ 
with T denoting time ordering. Expanding to the square 
order in the (i) , once the interaction is over the final 
state of the system is given by 
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where <&f = J 
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±mit <t>{x h t), and 6i = 
(Actually the last 



dt'mit 1 ) 

term in the expansion is of cubic order, because unlike 
the other cubic terms, it cannot simply be discarded at 



this stage.) When working in the computational ba- 

sis, {111, lit, IU, ITT, Til, TIT, TTI, TTT}, the detec- 
tors' nonnormalized reduced density matrix is given by 
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(3) 



Here we have employed the notation d""^ = 
(0| $f • • • $< |0), and C = (UU <0| 6, |0) ||||), with 
i, . . . , n = A, B, C and a, ...,( = ±. For simplicity, 
we have chosen temporally symmetric window- functions. 
Hence the amplitudes are all real. d^ + is the amplitude 
for a single photon emission by detector i, while dftjn 
is the amplitude for a single virtual photon exchange be- 
tween detectors i and j. The physical interpretation of 



the rest of the amplitudes should thus be clear. 

To prove that an iV-partite mixed state is truly ./V-fold 
entangled, it is enough to show that it can be distilled to 
a truly iV-fold entangled pure state 0, 0] . Having each 
of the detectors pass through a filter, which attenuates its 
"spin-down" component by a factor of r?, the detectors' 
nonnormalized reduced density matrix is in the compu- 
tational basis given by 
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Note that each of the components is written to its low- 



est nonvanishing order. The reason for this will shortly 
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become apparent. For L^ >> T, the overlap ampli- 
tudes, rf^iij) are negligible as compared to the emission, 

, and exchange amplitudes, dj^^. If we now take 
the window-function of detector C, and only detector C, 
to have a superoscillatory Fourier transform [l3L IT^ , of 
a form as in 0, E3, then by a suitable choice of the 
remaining two window-functions and the Lij the single- 
virtual photon exchange amplitudes involving detector 
C can be made arbitrarily larger than the rest pj, i.e., 

= d~£,\ ~S> all other amplitudes. In this limit, if we 
set rj 2 = = d~c\ and apply Wick's theorem [T^| . 

then it is readily seen that only amplitudes whose Wick 
decomposition contains terms consisting solely of single- 
virtual photon exchange amplitudes involving detector 
C survive the filtering. These amplitudes are of three 
types, representing single-virtual photon exchange pro- 
cesses d^, double emission processes di^ c , and the 
overlap between double emission processes by different 
detector pairs, df^^ (j ^ i). The detectors' reduced 
density matrix is thus given by 

/ 1 -1 -1 \ 
00000000 
00000000 

1 -1 0010100 ( r} 
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-100 1 1 00 

oooooooo 
Voooooooo/ 

This density matrix is pure and corresponds to the state 
-^jGIU) - IITT) - I TIT}), which, by means of local oper- 
ations on each of the detectors, can be transformed into 
a W state, ^ (|TU) + |1TI> + HIT))- The W state is 
truly tri-fold entangled. Since the detectors are initially 
nonentangled and remain causally disconnected through- 
out the interaction, this entanglement must have its ori- 
gin in corresponding vacuum entanglement, i.e., the vac- 
uum is truly tri-fold entangled. 

We now turn to the question of full nonlocality. The 
tri-partite W state violates the Svetlichny inequality > 
and is therefore fully nonlocal. Following the same rea- 
soning as above, we conclude that the correlations be- 
tween three arbitrarily separated regions of the vacuum 
do not admit a hybrid local-nonlocal hidden- variable de- 



scription. 

The generalization of the above analysis to any num- 
ber of distinct regions of the vacuum, N, is straightfor- 
ward. An TV-partite W state |3] can always be distilled. 
This can be seen as follows. Allowing for trivial modifi- 
cations resulting from the increase in the number of re- 
gions, we employ the same protocol as in the three-region 
case. We choose one and only one of the detectors, say 
detector A, to have the window- function with the su- 
peroscillating Fourier transform. Once the interaction is 
over, we have each of the detectors pass through a filter 
V 2 = d\jj = d++ = . . . = d^\ <N . (For any set of the 
L^, such a filter can be realized, as it is only the LiN that 
have to be taken into consideration, while the the ti (t) 
can always be suitably adjusted.) Only amplitudes whose 
Wick decomposition contains terms consisting solely of 
single-virtual photon exchange amplitudes, involving de- 
tector TV, survive the filtering. These amplitudes are of 
the same three types mentioned earlier. It is now only a 
matter of diligent book-keeping of the indices to convince 
oneself that, not excluding local operations, the filtering 
leaves us with an A-partite W state. It follows that the 
correlations between any A arbitrarily separated regions 
of the vacuum exhibit true A- fold entanglement. 

As regards nonlocality, the situation is different. For 
A > 4, numerical computations indicate that A-partite 
Svetlichny-type inequalities |Tj are not violated by the 
A-partite W state [ijj. This, of course, does not mean 
that the correlations between A > 4 distinct regions are 
not fully nonlocal. It may be that stronger inequalities, 
e.g. generalizations of the A-partite Svetlichny-type in- 
equalities incorporating more than two measurement set- 
tings per system |20j, will reveal the A-partite W state 
as fully nonlocal. Or a different detection scheme may 
make possible the distillation of states which are known 
to be fully nonlocal, e.g. multi-partite GHZ states. At 
present this question remains unresolved. 
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